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Initial states of glasma in high energy heavy ion collisions are longitudinal classical color electric 
and magnetic fields. Assuming finite color electric conductivity, we show that the color electric 
field decays by quark pair production with the life time of the order of QJ^, i.e. the inverse of 
the saturation momentum. Quarks and anti-quarks created in the pair production are immediately 
thermalized as far as their temperature is lower than Qs- Namely, a relaxation time of the 
quarks to be thermalized is much shorter than when j3~^ <C Qs- We also show that the quarks 
acquire longitudinal momentum of the order of Qs by the acceleration of the electric field. To discuss 
the quark pair production, we use chiral anomaly which has been shown to be very powerful tool in 
the presence of strong magnetic field. 
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I. INTRODUCTION 
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, \ Initial states of color gauge fields (glasma) produced immediately after high energy heavy- ion collisions have recently 
received much attention. The gauge fields are longitudinal classical color electric and magnetic fields; their field 
^SJ ' strengths are given by the square of the saturation momentum Qg. The presence of such classical gauge fields has 
J> i: been discussed on the basis of a fairly reliable effective theory of QCD at high energies, that is, a model of color glass 
00 I condensate (CGC)[l], Hj. It is expected that the decay of the glasma leads to thermalized quark gluon plasma (QGP). 

The glasma is homogeneous in the longitudinal direction and inhomogeneous in the transverse directions. Hence, 
we may view that it forms electric and magnetic flux tubes extending in the longitudinal direction. In the previous 
papersj^-Q we have shown a possibility that the famous Nielsen-Olesen instability [y| makes the color magnetic field B 
|n . decay. The possibility has been partially confirmed by the comparison between our results and numerical simulations 0, 
II]. Such decay is a first step toward the generation of QGP. 

On the other hand we have also discussed (ol-fTTj the decay of the color electric field; the decay is caused by 
Schwinger mechanism [l^ . that is, the pair production of quarks and anti-quarks. The mechanism has been ex- 
tensively explored [l 3] since the discovery of Klein paradox. Among them, the pair production in the expanding 
glasma has been discussed [l3l- A new feature in the glasma is that it is composed of not only electric field but also 
;J] \ magnetic field parallel to the electric field. Such a feature has also been explored. In particular, recently there are 
studies of the back reaction of the particles on the electric field under the presence of the magnetic field jl3, fisl [isj. 
The back reaction is essentially important for the decay of the electric field. 

Our originality [10| for the investigation of the decay is to use chiral anomaly. As is well known, the anomaly is 
effective when collinear magnetic and electric fields are present. This is the very situation in the glasma. When we 
use the chiral anomaly, we can discuss Schwinger mechanism without detail calculations [TsUTsj of wave functions but 
simply by solving classical anomaly equation and Maxwell equations. In particular, when the strong magnetic field is 
present, the anomaly is much simplified because the quarks are effectively massless and only relevant states are ones 
in the lowest Landau level. ( Both \/ gE ~ Qs and ^JgB ^ Qs in the glasma are much larger than mass of quarks. 
) Since the motions of the quarks in transverse directions are frozen, only possible motion is along the longitudinal 
direction. Thus, the anomaly equation is reduced to the one in two dimensional space-time. With the simplification, 
we can find important quantities even in complicated realistic situations for which the investigations have not yet 
performed. Actually, we have discussed the decay of axial symmetric electric flux tube by taking account of the 
azimuthal magnetic field around the tube. The field is generated by the current carried by the pair created quarks 
and anti-quarks. Although the electric field loses its energy by the pair creation and the generation of the azimuthal 
magnetic field, it never decays. The field oscillating with time propagates to the infinity in the axial symmetric 
way[l3- This is because the quarks are free particles and there is no energy dissipation. ( In the case of homogeneous 
electric field, the field simply oscillates with time. ) 

In this paper we wish to discuss the decay of the color electric field E by taking account of energy dissipation in 
heat baths. The dissipation arises due to the presence of finite electric conductivity. Namely, the pair production 
generates color electric current J, which dissipates its energy owing to the interaction between the quarks and the 
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surrounding; the surrounding is composed of quarks and gluons. Actually, we assume Ohm law J = aE with electric 
conductivity a. The conductivity is calculated by using Boltzmann equation in the relaxation time approximation. 
In the approximation a relaxation time is obtained by calculating electron's scattering rates. Then, we can show 
that the quarks are thernialized immediately after their production as far as their temperature is much smaller than 
Qs ~ V eE ; the relaxation time of a slightly deformed momentum distribution of quarks becoming the equilibrium 
Fermi distribution is much shorter than the life time ^ Qj^ of the field. As numerical calculations have shownfl^, the 
longitudinal momentum distribution of the free particles produced in the vacuum is almost equal to the equilibrium 
one, that is Fermi distribution at zero temperature. Thus, even in non zero but low temperature, the momentum 
distribution is nearly equal to the equilibrium one. Our relaxation time approximation in Boltzmann equation may 
holds in such a situation. Therefore, owing to the energy dissipation by the scattering between electrons and positrons, 
the electric field decays and never oscillates. 

For simplicity, we examine homogeneous electric and magnetic fields of U(l) gauge theory instead of SU(3) gauge 
theory. Thus, we use terminology of electron or positron instead of quark or anti-quark. The generalization of our 
results to the case of SU(3) gauge theory is straightforward done simply by assuming maximal Abelian group of 
SU(3)[l^. We assume that both the electric and magnetic fields are much larger than the square of the electron mass. 
Thus, they are taken to be massless. 

In the next section we explain how the chiral anomaly is useful for the discussion of Schwinger mechanism. We 
apply the anomaly to the discussion of the pair production with energy dissipation in the section (3). In the discussion 
we use the conductivity of electron-positron gas under strong magnetic field. The detail calculation of it is shown in 
the appendix. 

II. CHIRAL ANOMALY AND PAIR PRODUCTION 

We first make a brief review of the utility of the chiral anomaly. As we will show below, only by simple calculations 
we can obtain most of important results of Schwinger mechanism, which was derived previously by detail calculations. 

Under the homogenous magnetic field B, massless fermions with charge e > occupy the states characterized by 
Landau level and momentum p parallel to the magnetic field. Their energies are given by 

En = \/p'^ + 2NeB (parallel) and En = ^^p'^ + eB{2N + 1) (antiparallel), (1) 

where integer N denotes Landau level. The term of "parallel" ("anti-parallel") implies magnetic moment parallel 
(anti-parallel) to B. The magnetic moment of electrons (positrons) is antiparallel (parallel) to their spin. Thus, 
electrons (positrons) with spin anti-parallel (parallel) to B can have zero energy states in the lowest Landau level; 
their energy spectrum is given by En=o = \p\ > 0. On the other hand, the other states cannot be zero energy states. 
Hence only relevant states for the pair production are states with energy En=o — \p\ in the limit of strong magnetic 
field. Hereafter, we consider only such states. 

The equation of the chiral anomaly for the states leads to 

dtinR~nL)^^m-B (2) 

where { ) denotes number density of right ( left ) handed chiral fermions; tij^^l = (^7o(l ± 75)^')/2 in which 
the expectation value is taken by using the states in the lowest Landau level. We have used spatial uniformness of 
the chiral current J5, that is, divjs = 0. 

We note that electrons move to the direction anti-parallel to E while positrons move to the direction parallel 
to E after their pair production. Therefore, both the electrons and positrons have right handed helicity when E 
being parallel to B. ( Right or left handed helicity means the state with momentum parallel or antiparallel to spin, 
respectively.) Since chirality is equivalent to helicity of massless fermions, n/j = 2n and = where n is the number 
density of electrons. ( When E is anti-parallel to B, nji = and ul = In. ) Therefore, the anomaly equation 
describes the pair production rate of electron and positrons under E and B, 

2 2 
dtn=-^E(t)B for E{t) > and dtn^^-^E{t)B for E(t) < (3) 

For example, we suppose that the electric field E{t) > is switched on at i = 0. Then, solving this equation with 
initial condition n(t = 0) = 0, we can obtain the number density of electrons fisl; n{t) = ■^E{t')Bdt' . When E is 
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constant, n increases linearly with time. In this estimation, the back reaction of the produced particles to the electric 
field is not taken into account. 

When we take account of the back reaction, we can find that the electric field oscillates with time. The back reaction 
is that the energy of electric field decreases owing to the acceleration of the charged particles. It can be included by 
solving a Maxwell equation, dtE = — J along with the anomaly equation. To solve the equations, we need to express 
the electric current J in terms of n. This can be done by considering the energy conservation. Namely, the total 
energy of the electric field and the charged particles e is conserved, 

dt + ^^"^^O " + dtn{t)pF + n{t)eE ^ E{^J + 2en{t)) = 0, (4) 

where the energy density of electrons or positrons is given such that t{t) = n(t)pp{t)/2 with Fermi momentum 
PF{t) — Jq eE{t')dt' — n(t)4:n'^ / eB for E > 0. Namely, we assume the energy of free electrons and positrons. Their 
momentum distributions are supposed to be no oc d{pF{t)—p)0{p) for positrons or no oc 9{pF(t)+p)9{~p) for electrons. 
We have used the formula dtn{t)pF — dtn{t)n{t)4:TT'^ / eB = en{t)E in eqQ. Thus, from the equation we obtain the 
electric current J = 2en for E > 0. Obviously, there is no energy dissipation. 

Using the Maxwell equation and the anomaly equation we derive the equation of the number density, 

dMt) + ^n{t)^0. (5) 

We can see that after the switch on of E, E gr adually decreases with time, while n increases. Eventually, these 
quantities oscillate with the frequency e-\/ eB jl-K^ ( plasma oscillation ) as shown in Fig.l. The oscillation can be 
understood in the following. 
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FIG. 1: number density n (dash) and electric field (line) with arbitrary scale 

When i? ( > ) parallel to B is switched on at i = with n(t = 0) = 0, the number n of the particles increases while 
E decreases. When E vanishes, n takes the maximum. After that E changes its sign, namely, E becomes negative ( 
anti-parallel to B.) Then, the particles are forced to move to the direction inverse to the one in the case of E' > 0. 
Thus, the pair annihilation occurs and n begins to decrease. The energy of the particles is transmitted to the energy 
of the electric field. Actually, according to the anomaly equation, \E\ becomes large with time. When n vanishes, 
\E\ takes the maximum ( E = —E{t = 0) ). Subsequently, the pair creation occurs and n increases according to the 
anomaly equation(|3]) but with E < 0. This is a plasma oscillation shown in the figure. The energy of the electric 
field is transmitted to the particles and then the energy of the particles goes back to the electric field. The plasma 
oscillation is a natural consequence of the absence of the dissipation. 

We should make a comment on the scales involved in the discussion. The above result ( plasma oscillation ) holds 
in the limit of eB ^ eE where only the states in the lowest Landau level are relevant. On the other hand when 
eB eE, states in higher Landau levels become important since the energy |p_F(t)| = | Jq dt'eE\ of electrons can be 
much larger than V eB. In the glasma we have both eB and eE ~ ^ ™cicctron ( square of electron mass ) . But 
states in higher Landau levels are still not important according to the following reason. When some of electrons in 

the lowest Landau level are accelerated and their energies reach the energy |pf| = /g ' dt'eE ~ eEQ^^ ^ VgB, they 
may make a transition to a state in a higher Landau level with energy = VgB. Such electrons are the ones produced 
in the very early stage of the pair production. The number of the electrons is still a fraction. Most of electrons whose 
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energies are less than \J eB, occupy the lowest Landau level. It has recently been shown[15j that the decay of the 
electric field is mainly caused by particles in the lowest Landau level. Thus, our result is approximately correct. 

In this way, we obtain the number density of electrons simply by using the chiral anomaly. It is important to note 
that the anomaly equation involves all of quantum effects. Thus, by simply solving the classical anomaly equation, 
we can obtain the quantities resulting from the purely quantum effect, i.e. Schwinger mechanism. In the standard 
manipulation [isl] of Schwinger mechanism, we need to find wave functions of electrons under the electric and magnetic 
fields in order to obtain such quantities. When the time dependence of the electric field is complicated, it is very 
difficult to obtain the wave functions. On the other hand, the anomaly equation is a very simple and powerful tool 
for investigating Schwinger mechanism. 



III. PAIR PRODUCTION WITH ENERGY DISSIPATION 



Up to now, we have not taken into account energy dissipation of electrons and positrons produced by the electric 
field. We proceed to analyse the Schwinger mechanism in heat bath, especially, the case with energy dissipation. 
We assume that the particles are immediately thermalized after their production. ( Later we see that relaxation 
time, within which slightly deformed particle distribution becomes equilibrium distribution, is much shorter than the 
typical scale Qj^ ~ ^/TJeE as far as the system is in low temperature. As numerical calculations have shown [T3|, 
the momentum distribution of the produced particles is nearly equal to the Fermi distribution at zero temperature. 
) Thus, the number density n is given by 



+00 

^^exp(/3(t)(|p|-/i(t))) + l' 

with ao = eB / 27r^ , where /3 = 1 /T is the inverse of temperature T and /i is chemical potential of electrons or positrons. 

The number density n depends on two parameters, /3 and fi, which may be functions of time. Hereafter, we 
consider only the case that the chemical potential /x vanishes due to non-conservation of electron number. It is caused 
by frequent pair annihilation or creation of electrons and positrons after their production. Since initial states have 
no lepton number, the chemical potential vanishes. ( This corresponds to the situation in the glasma where the 
pair annihilation or creation of quarks and anti-quarks frequently take place after the pair creation by Schwinger 
mechanism. ) Thus, the number density is only a function of the temperature and is given by n = ao/3(i)~^ log 2. 
Similarly, the energy density is given by 



n{t) = ao 



e{t) = ao / dp J = coaoPit)-^ (7) 

Jo exp {I3{t)p) + 1 

with Co = xdx/{exp{x) + 1) = 7r^/12. Since the temperature f3~^ represents an average energy of the particles, it 
is reasonable that we have e ~ the energy density is equal to the number density times the average energy of a 

particle. 

The electric field produces electrons and positrons and generates their electric currents. When the electric current 
flows, the current dissipates its energy due to the finite electric conductivity of these particles. Thus, the electric 
field dissipates its energy. Here we may assume that the electric current J satisfies the Ohm law J = cr (/?)£' 
where a denotes electric conductivity. Obviously, the Ohm law represents the energy dissipation of the electric field 
dt{E^/2) = EdtE = -EJ = -aE^ < 0. 

As we will show in the appendix, the electric conductivity a of the charged particles in the lowest Landau level is 
given by 



where /o = l/(exp(|p|/3) -I- l) denotes the momentum distribution of the charged particles. The conductivity depends 
on the relaxation time t(p); the time needed for a momentum distribution to become the equilibrium one fo{p)- The 
relaxation is realized by scattering of the charged particles so that it can be expressed by their scattering amplitudes. 
The detail calculation shows that r {eBl3)~^ in the limit of low temperature (3~^ <^ VgB. ( This fact makes valid 
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our assumption of the immediate thermalization in the low temperature where the relaxation time is small. ) Thus, 
the conductivity depends on the temperature such as (j{f3) ~ bo/3~^{t) in the limit ( 60 is a numerical constant ). 

Thus, using the Maxwell equation dtE — —J — —aE and the anomaly equation dtn = ( 7 = e'^B/'in'^ ), we 
find that 

d?\ogE{t) = -dta = -bodj^^) = -bJ^^) 

V flQ log 2 / \ao log 2 / 

Setting X = \og[E{t)/E{t — 0)), we obtain the equation of the electric field 

d?X{t) + ^°y^('^^= °) exp(X(0) ^ (10) 

with initial conditions X{t = 0) = and dtX{t = 0) = 0. The number density n{t) can be obtained by integrating 
the solution E{t) = E{t = 0) exp(X(t)) such that n{t) = 7 J* dt'E[t'). 

The equation (flO)) describes a motion of a classical particle with the potential energy exp(X); X denotes the 
coordinate of the particle. Starting at X = 0, the particle goes down to —00 along the slope of the potential. It means 
that E cx exp(X) — )■ 0. Namely, the electric field vanishes with the life time of the typical scale 1/ \/e E{t — 0) ^ l/Qs- 
( The life time of the electric field is sufficiently short to be consistent with phenomenology [iTj of QGP. ) Since 
n{t) =7/0 dt'E{t') = aoT(t)log2, the number density n or the temperature T monotonically increases with time 
and reaches its maximum. This is contrasted with the case of no dissipation in which both E and n oscillate. These 
behaviors are shown in the Fig. 2. 
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FIG. 2: number density n (dash) and electric field (line) with arbitrary scale 



We can roughly estimate the attainable number density n, temperature T and average energy of the particles ^ e/n 
in the pair creation. 



e'B r e^B E{t ^ 0) 



47r2 7o 47r2 ^eE{t = 0) 



y^eE{t = 0) and - - T (11) 



n 



1 which corresponds to gauge coupling constant f;^ ~ 1 in QCD with the energies at RICH or LHC. 
We should make a comment. The result has been derived in the limit of the low temperature. The relaxation time 
T oc /3~^ is sufficiently short in the limit so that we can use the equilibrium momentum distribution. The temperature 
increases with time from zero after the switch-on of the electric field. Therefore, the result is reliable in the early stage 
of the pair creation. However, when it becomes comparable to \/gB, the assumption does not hold. Then, we cannot 
use the formula a ~ t oc f3~^ and the equilibrium momentum distribution. Hence, the result of the saturation 
in the number density at the late stage of the evolution cannot be reliable. But, as the temperature increases and 
becomes comparable to y/gB, the conductivity is almost constant in time. Then, we find that E — Eq exp(— crt) from 
the equation dtE — —aE and that n = tiq — jcjEq exp(— cri) from the anomaly equation. The number density is 
saturated and never oscillate. ( To obtain this result we do not use the equilibrium momentum distribution. Namely, 
the assumption of the particles being immediately thermalized is not necessary. ) Actually, when the temperature 
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or the average energy of the particles is comparable with yfgB, the electric conductivity may depends weakly on the 

time because only dimensional quantity in a is \/gB or the average energy. Consequently, we find that the number 
density ( temperature ) gradually increases in the early stage and is saturated at the late stage. It never oscillates 
unlike the case with no energy dissipation. 

To put it another way, in the early stage when the pair creation starts, the particles are thermalized soon because 
the relaxation time is short. Thus, the temperature can be defined. The temperature or the average energy of the 
particles is much less than eE ~ Qs- But in the stage when relaxation time is not sufficiently short for the particles 
to be thermalized soon, the temperature cannot be defined. In the stage all of the typical energy scales are of the 
order of Qs. Thus, the electric conductivity is also of the order of Qs- This leads to the above result of the saturation 
n = hq — ^aEo exp(— crt). The stage is not early stage but may be late stage in the pair creation. In this way we 
understand that the results shown in Fig. 2 are reliable. 



We have shown that the color electric field decays into the quarks and anti-quarks under the color magnetic field. 
They occupy only the states in the lowest Landau level. We have also discussed that the average energy of the 
particles is of the order of the saturation momentum Qs- Namely, the Fermi momentum, which is the longitudinal 
momentum, is given by Qs- The particles obtain the longitudinal momentum with the acceleration by the electric 
field. We remind you that the glasma has only transverse momentum of the order of Qs and has no longitudinal 
momentum owing to the homogeneity in the direction. Thus, initially there is anisotropy in momentum distribution. 
But owing to the decay of the electric field, the particles have longitudinal momenta of the order of Qg- Thus, the 
longitudinal and transverse momenta become almost identical to each other. This is a step toward the isotropy of the 
momenta or thermalization of the quarks. 

The life time of the electric field has been found to be of the order of Q~^, while that of the magnetic field by a 
decay mechanism of Nielsen Olesen instability is much longerQ than Qj^. We may speculate that the magnetic field 
also decay rapidly due to Schwinger mechanism. This is because azimuthal electric field is generated by the expansion 
of the longitudinal magnetic flux tube and this electric field decays by Schwinger mechanism. The expansion in 
the transverse direction of the magnetic flux tube is simply a result of Maxwell equations 0, 0]. Similarly, azimuthal 
magnetic field is also generated by the expansion of the longitudinal electric flux tube. Thus, both azimuthal magnetic 
and electric fields are generated around the longitudinal flux tubes. Then, the azimuthal electric field decays in the 
way as we have shown. The field strengths of these longitudinal and transverse ( azimuthal ) fields have been shown 
to become the same order of magnitudes at about the time Qj^ after the high energy heavy ion collisions. In this 
way, the longitudinal magnetic flux may decay rapidly through the decay of the azimuthal electric field generated by 
the expansion of the magnetic flux. Such a rapid decay of the magnetic field as well as the electric field is favorable 
for establishing early thermalization , 1 7j of QGP. 

To summarize, we have shown how the electric field decays by Schwinger mechanism in the presence of the energy 
dissipation. Namely, the electric current carried by electrons and positrons dissipate its energy owing to the Ohm law, 
which results in the increase of the temperature with time. Our result holds only in the limit of low temperature, in 
which the relaxation time has been shown to be much short; the thermalization of produced particles can be realized 
immediately after their production. However, we have also discussed that our result is approximately correct even 
when the temperature or the average energy of the particles is of the order of Qs- 



We show that the electric conductivity behaves as cr cx /3 ^ in the limit of the low temperature — oo. We 

calculate a in the relaxation time approximation of Boltzmann equation. 



IV. DISCUSSION AND SUMMARY 



Appendix A: electric conductivity (J 



PidpJ = PzdpJ = eEdpf 



f-fo 



(AI) 



T 



with p = Pz, where we have used classical equation of motion, p = eE + ep/\p\ x B. t denotes relaxation time and 
/o does Fermi momentum distribution with the vanishing chemical potential; /o = l/(exp(|p|/3) + l) . / represents a 
deformed momentum distribution caused by the effect of the electric field. Since the electric current J is given by 



7 



/-f-oo /"f-oo p-i-oo p-roo 



we obtain the electric conductivity 



a 



+00 

dpT{p)h{p){l- fo{p)). (A3) 



47r2 

The relaxation time r(p) is given by the inverse of the scattering rate W{p) between electrons and positrons, 

e~{p) + er[q) e"(fc) + er [1) and e"(p) + e+(g) -> e"(fc) + e+(/) (A4) 

1 1 1 1 

t{p) = + + — + 



W^dir(e-e-) Wex(e-e-) l^crossle-e-) T^an(e+e-) 

111 
^Tydir(e+e-) ^ Wcx(e+e-) ^ Tycross(e+e-) ^ 

where we denote longitudinal momentum of each particle ( parallel to i? ) as p, fc,and I. We use the notations 
like Wdir, etc. to distinguish the contributions of different Feymann diagrams, which was used in the paper [18J. We 
also use the notations VF(e"e~) and W{e^e^) to distinguish the contributions of the processes e^e^ -> e~e~ and 
e+e" — >■ e+e". 

Our concern is scattering of the massless fermions occupying only the states in the lowest Landau level. Their wave 
functions are given by 

/eB\i/4 ^ _ . . ^ f eB 



with Ep — zL\p\ and Dirac spinor u{p), 



exp{-iEpt + ip2y + ipz) exp(^-^{x - (^6) 




where we use the gauge potential A — (0, xB, 0) of the magnetic field B — (0, 0, B). Ly and denote spatial lengths 
in the direction of y and z axis, respectively. The states are characterized by two momenta p and p2. Using the 
formulae, the scattering amplitudes can be calculated in the following way. For example, the scattering amplitude 11 
of e~(p) + e~{q) -> e~(fc) + e~(Z) associated with Vydir(e~e~) is obtained by estimating the formula, 

d^xd^yJ^ix;p,k)J'^iy;qJ)'-^^^^^^^^^^^^^ (A8) 

with Ji^i_{x;p,k) = e^p(x)7^t^fc(a;) and q'x — q'^x'^ — q[x^ — q2X^ — Qa^^ ■ Here we also used another notations 
t = x'^ ,x — x^ ,y ~ x'^ , z = x'^ in ea (jA6p . Using the amplitude 11, we define the scattering rate Wdir(e~e~), 

„ o LtLtdqydqdkydkdlydl 

dWdn-ie-e-) =. \Ufi2TTfS{po + go - fco - Wpy + qy - ky - ly)5{p + q-k^ I) ^ J ^. (A9) 

[Ztt 



The other scattering rates can be obtained in the similar way. Consequently, we find the formulae of W{p) which are 
fairly simplified since the fermions are massless. 



8 



g4 r + oo 

Wdir{e^e^;p) = —— / dqdkdldq2dk2dl2S{p + q- k- l)6{p2 + q2- k2- l2)5{\p\ + \q\ - \k\ 



-OO 

1 Z^™-^ ^'^^^'^^'^ /•^^.«^p(-^-ot)A^ 



X ((|p||fc| +pk){\q\\l\ + ql) - {\k\p+\p\k){\l\q + l\q\)f fo{q){l - /o(fc))(l - /o(0). (AlO) 

where Z2.^ = fc2 + 2(|p||fc|-pfc), 



g4 /•+00 

Wecc(e~e~;p) = — — — / dqdkdldq2dk2dl2S{p + q - k - l)S{p2 + q2 - k2 - l2)S{\p\ + \q\ - \k\ - \l\) 



2 



X ((H|;| +pO(kPI + qk) - i\i\p + \p\i)i\k\q + k\q\)f /o(g)(i - /o(fc))(i - /o(0), (Aii) 

where Zl = q + 2{\p\\l\- pi), 

g4 r+oc 

Wcross{e~e~;p) = / dqdkdldq2dk2dl2S{p + q-k- l)S{p2 + 92 - - ^2)^(1^1 + \q\ - \k\ - 



(blkl|fc||/|) 



2 



X ((H|fc| + (|fc|p + Hfc)(|;|(7 + ;|(7|)) 

X ((H|Z|+pO(l9l|fc|+3fc)-(l^|P+blO(|A:k + fc|g|)) /o(g)(l-/o(A:))(l-/o(0), (A12) 



and 



g4 /-i-oo 

Wa„(e+e~;p) = -— j / dqdkdldq2dk2dl2S{p + q-k- l)5{p2 + q2-k2- l2)S{\p\ + \q\ - \k\ - \l\) 
1 .expf-(^±i^^^)f,/^exp(-^)- 



(b||g||fc||/|)2Z2„ -V eS nV2 eS Wo x/2^ / 

X (blk|-OT)'(|fc||/|-fcOVo(5)(l-/o(fc))(l-/o(0), (A13) 

where Zl^ = {k2 + h? - 2{\p\\q\ - pq), 



r+oo 

4(27r)i" 



Wdir{e+e ;p) = ,/o_^4 / dg(ifc(iWg2C«fc2C^^2'^(p + 9 - - 0<5(P2 + (72 - ^2 - ^2)<5(bl + |g| - |fc| 



{{\p\\k\ +pk){\q\\l\ + gO - {\k\p+\p\k){\l\q + l\q\)f /o(g)(l - /o(A:))(l - /o(0), (AM) 



g4 (-+00 

Wex(e+e~;j3) = / dg'(ifcdWg'2*2rfZ2^(p + g - A; - 0^(P2 + 92 - ^2 - /2)5(|p| + |g| - |fc| 
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i\pMk\WZi 

{Mi\+pi)Mk 



exp 



qk) 



^2 ~ ^ex 

eB 



(A/fexp(- 



Zex\k2 + hi - 

eB 



dx 



exp(— a; — 



\vmk\q + k\q\)f fo{q)(l - /o(fc))(l - /o(0), 



(A15) 



Wcross{e~^e ;p) = / dqdkdldq2dk2dl2S{p q - k - l)S{p2 + g2 - - l2)S{\p\ + |g| - \k\ - \l\) 



1 .„p(_M±t£l.^)(,/Iexp(-5^lttl±M,_ '^^^SS?-) 



V 2eB yVV2 eB ' Jo ^ 



X 



((H|fc| +pfc)(|g||Z| +qO - {\k\p+\p\k){\l\q + l\q\)) 



X ((H|/|+pO(kl|fc| + 'zfc)-(l'b + blO(|fck + fckl)) /o(<z)(i - /o(fc))(i - /o(0), (ai6) 

where the deha functions represent momentum and energy conservation. Ah of the formulae involve the integration 
variables, g, k, I, q2, ^2 and h- 

Now we explain that the electric conductivity behaves as cr oc fi"^ in the limit of gB0^ — ?> oo. In the integrals of W , 
we rescale all momenta such as e.g. p = [3~^p. Then, we find W{p) = [3~^W{p) in which dimensionless quantity W{p) 
depends on (3 and eB through the dimensionless quantity eB(3^. Using the delta functions we can trivially perform 
the integration of the variables g, q2 and I. Among remaining integrations of fc, ^2 and h, the integration of k is finite 
at |fc| = oo due to the factor fo{q = —p + fc + F) ~ exp(— for |fc| — >• oo. On the other hand, the integrations of ^2 
and h are finite at \k2\ — oo or Ih] = oo owing to the factors exp(— fc|/ei?/3^) or exp{—ll/eBl3'^) and l/^jir ^/M 
or l/T^. Hence, W behaves such that W{p) = f3-^W{p) oc (3-'^eB(3'^ in the hmit of eB/S'^ oo. Therefore, 

<T oc PeB J dprlp) oc PeB J dpl/W{p) oc 

As we have stated before, the relaxation time r ~ {eBf3)~^ ^ {T/Qs)^/Qs is quite shorter than the typical time 
scale l/^yeE{t = 0) ^ Qj^ in the period of T/Qg ^ 1. Thus, we may think that the particles produced by Schwinger 
mechanism are immediately thermalized at least in the early stage of the electric field decay. 
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